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We investigate the severity of the sign problem in a random matrix model for 
QCD at finite temperature T and baryon chemical potential [i. We obtain analytic 
expression for the average phase factor - the measure of the severity of the sign 
problem at arbitrary T and (i. We observe that the sign problem becomes less severe 
as the temperature is increased. We also hnd the domain where the sign problem is 
maximal - the average phase factor is zero, which is related to the pion condensation 
phase in the QCD with finite isospin chemical potential. We find that, in the matrix 
model we studied, the critical point is located inside the domain of the maximal 
sign problem, making the point inaccessible to conventional reweighting techniques. 
We observe and describe the scaling behavior of the size and shape of the pion 
condensation near the chiral limit. 



I. INTRODUCTION 



The phase diagram of QCD at finite temperature and baryon density has been a subject 
of intense interest during recent years (see, e.g., [1] for review). On the one hand, the 
experimental heavy-ion collision programs, such as those at CERN SPS and RHIC, as well 
as planned FAIR facility, demand reliable quantitative understanding of the phase diagram 
of strongly interacting matter created in those collisions. On the other hand, deriving the 
relevant predictions from the first principles of QCD is a formidable challenge, because the 
phenomena associated with phase transitions occur in the domain where the QCD coupling 
is not small. 

One of the features of the QCD phase diagram of particular interest to heavy-ion collision 
experiments is the QCD critical point. The existence of such a point - an ending point of 
the first order chiral transition in QCD - was suggested a long time ago , and the 
oroperties were studied using universality arguments and model calculations more recently 



(see Ref. 



6| for review). The experimental search for the critical point using heavy ion 



collisions has been proposed in [7] . It is apparent that theoretical knowledge of the location 



2 



of the critical point on the phase diagram is important for the success of the experimental 
search. 

The time-tested approach to non-perturbative problems in QCD is the numerical lattice 
Monte-Carlo calculations. This approach, very powerful at studying QCD thermodynamics 
at zero baryon density, runs into the sign-problem at nonzero baryon density. The lattice 
calculations are based on reinterpreting the QCD partition function as a partition function 
of a classical statistical system, with energy given by the Euclidean action of QCD. This 
action involves the logarithm of the fermion determinant, which is complex for any nonzero 
value of the baryon chemical potential \lb- The Monte Carlo importance sampling technique 
uses the exponent of the action as the measure of importance and fails because the action 
is complex. 

Several approaches to QCD at finite baryon density are being developed, using various 
techniques to circumvent or tame the effect of the complexity of the fermion determinant 
and locate the QCD critical point Q, <J Q U, 12] (see also Refs. [3, 14, Jjj, Q for 



reviews). As the sign problem worsens with increasing the hope is that the QCD critical 
point is located at sufficiently small \lb, where the sign problem can still be controlled. In 
many cases it is difficult to judge reliably, either a priori or a posteriori, what the range of 
validity of the results are, in terms of how large /ig can be before the sign problem is out of 
control. It is therefore necessary to understand better the severity of the sign problem, and 
its dependence on the variables such as temperature T, baryon chemical potential fiB and 
quark mass m. 

In this paper we use a random matrix model of QCD to assess the severity of the sign 
problem as a function of T, and m. A similar study at T = has been reported in 



Ref. 



13, 



18| . Here we shall present analytical 1 results for the random matrix model at 
nonzero temperature - the regime most relevant for the heavy-ion collision experiments and 
the lattice studies aimed at discovering the QCD critical point. 

As a quantitative measure of the severity of the sign problem we consider the complex 
phase e ld of the fermion determinant det D, averaged over gauge field configurations of the 
phase-quenched theory: 



1 A numerical study at nonzero T has been reported in Ref. 19 1. 
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where D is the Euclidean space Dirac operator in a given gauge configuration: 



B = 'f(d li -A lt )+m + iiy . 



(2) 



and /i is the quark chemical potential: /i = /zb/3. 

The average in Eq. ffll) . denoted by (. . is taken over the gauge field configuration 

ensemble with the phase of the determinant removed (quenched), making the measure of 
path integration manifestly positive: e~ SYM | det B| 2 . This phase-quenched theory can be 
viewed as a theory with 1 quark and 1 conjugate quark or, due to (detB(/z))* = detO(— /i), 
two quarks with opposite chemical potentials, i.e., QCD at finite isospin chemical potential 
fi! = 2/x (see, e.g, 



20 



2l|) 



The average phase factor R can be recast as the ratio of two partition functions 

<(det©) 2 )o Zi+i 



17 



19|: 



R = {e 2t0 ) 



1+1* 



(3) 



(|detD|2) Z 1+1 .' 

where (. . .) denotes average over gauge configurations in a theory without quarks (with 
pure Yang- Mills measure e~ SYM ). 



II. RANDOM MATRIX THEORY (RMT) AND QCD PARTITION FUNCTION 



The random matrix model 



The Chiral Random Matrix Theory 22| approximates the QCD partition function by an 



integral over random matrix ensemble. We introduce temperature as first done in Ref. [23] 
and chemical potential as in Ref. 24|. The resulting random matrix model has been used in 
Ref. [5( to study the QCD phase diagram. The partition function in the model is given by: 



Zn, 



VXe 



det^D = (det^B)* 



(4) 



where B is the 2N x 2N matrix approximating the Dirac operator: 

( m iX + C^ 

B = 



m iX + C 
z'JSft + C m 



(5) 



with 



C = fil N + iT 



( In/2 ^ 



-1 



(6) 



N/2 J 
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X is an JV x JV complex random matrix; Ijv is the JV x JV identity matrix. The deterministic 
matrix C defined by Eq. ([6]) accounts for the effect of the chemical potential \x and of the 
2 smallest Matsubara frequencies: +irT and — 7rT. For simplicity, we absorb the coefficient 
7r into T. The integration in Eq. (J3J) is over the real and imaginary components of the 
matrix X: VX = Y] N i dX u dXh . 

All quantities appearing in the random matrix model are dimensionless, which is achieved 
by using the appropriate units as discussed in Ref. The choice of the dimensionful units 
will not be consequential for our study. 

The Dirac determinant can be written as a Grassmann integral: 



det^D = / IJ ^ f R T>^{ ex P 



E 



m 



iX + C 



\ 



V 



iXi + C 



m 



I 



\ *i ) J 



(7) 



where the Grassmann integration is over the spinors: T>^ = Yli=i dip?* . 

The integration over random matrix X is Gaussian and leads to a 4-fermion interaction: 
(ipi* r 4 )9 L i ipR. ipRj J here i and j indicate dimension of random matrix X while / and g indicate 
the flavors. 

Following the logic of the Hubbard- Stratonovich transformation, 4-fermion interaction 
can be rewritten using fermion bilinears with the help of a new auxiliary iVf x JVf complex 
matrix A (flavor matrix). Performing the Grassmann integration one then obtains 



VAe (-NtrAAl) d ^ 



^ A + m fi + iT^ 
yfi + iT A* + m j 



det^ [T -> -T 



(8) 



We shall specialize to JVf = 2 quark flavors. The integral in Eq. (jHJ) is performed over 
2 x JV f x JV f = 8 variables which are real and imaginary parts of the elements of the complex 
JVf x JVf = 2 x 2 flavor matrix A. We shall define potential Qi + i(A) as 



Z 1+1 = J VAe~ Nn ^ A \ 



(9) 



i.e., 



n 1+1 (A) = Tr 



AA^ --\n{[{A + m){A^ + m)-{fi + iT) 2 ] x [T^-T]} 



(10) 
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B. Phase quenched partition function 

The phase-quenched partition function is given by 



'i+i* 



J VX e~ N trXXt detD detD* = ( |det©| 2 ) x ■ (11) 



Following the steps outlined above, the phase quenched partition function Zi + i* can be 
written as 

/ 



Z 1+1 * = J VAe(~ NtTAA ^det 



det^ [T^-T] (12) 



A + m fir 3 + iT 

, fJLT 3 + iT A^ + 171 

= I VAe- Nn ^* (A \ (13) 



Here t% is the Pauli matrix, and we defined another potential f2 1+1 *(A) — a function of a 
2x2 complex flavor matrix A. For a generic matrix A, fi 1+1 .(A) ^ Qi + i(A), due to the 
presence of the Pauli matrix in Eq. (|12p . 



C. Thermodynamic limit N — > oo and the solution of the model 

In the thermodynamic limit iV — > oo, the random matrix partition function Zi+i can be 
calculated analytically using the saddle-point approximation by minimizing f2 1+1 (v4) with 
respect to A. The minimum is given by a multiple of the unit matrix: A — at with a - real. 

For the 1 + 1* theory the minimum of Qi + i*(A) is also given by a multiple of the unit 
matrix, except for a region on the phase diagram where another, deeper minimum is given 
by a non-diagonal matrix A 2^]. This breaks the U(l) (ts isospin) symmetry of the theory 
in this region and is associated with the pion condensation in QCD. 

Outside of the pion condensation region, i.e., when the minimum of is given by a 

multiple of the unit matrix, the minimum values of the two potentials coincide. Indeed, for 
any real a, f2i + i(al) = f2i + i*(al) = Q(a), where we defined 

fi(a) = 2a 2 - ln{[(a + mf - {iT + ^) 2 ][(a + mf - (iT - /i) 2 ]} . (14) 

Then the saddle-point value of A — at for both and Zx+i* is determined by minimizing 
the potential Q(a) with respect to a: 

( a + m )[( a + m )2_ / ,2 +T 2 ] 

[(a + mf - /x 2 + T 2 ] 2 + VT 2 ' 1 ' 
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III. AVERAGE PHASE FACTOR IN PHASE QUENCHED THEORY 



A. General result: arbitrary m, fi and T 



Since the leading exponential behavior of the partition functions -Zi+i and -Zi+i* is the 
same (see previous section), it will cancel in the ratio R. Therefore, we have to take into 
account the preexponential factors, which are determined by the second order derivatives of 
the potential function Qi + i(A) and (A) with respect to all elements of flavor matrix A: 

J Q 



-nq q (A) 



(16) 



A=A S 



where Q indicates the quark content of the theory, 1 + 1 or 1 + 1*, and 

det £l n = det ' 



dA a dA p J ' ^ 
The indices a and (3 run through eight values labeling eight independent components of the 
complex 2x2 matrix A: A a and Ap = (A n , A\ x , Ai 2 , A* 2 , A 2 i, A* 21 , A 2 2, A^). Evaluating 
determinants in Eq. (TlTl) at the saddle point A = A sp we find, using notations given below 
in Eq. (|20|): 



det — [x — y ] and det — [x — y ] [u —v 



,2l2 



(18) 



Since Qi + i(A 8p ) = Q\ + i*(A sp ) outside of the region of pion condensation in 1 + 1* theory, 
the exponential factors in Eq. ffl6|) cancel in the ratio R and the average phase factor is 
given by 



R = (e 2 ^) 



Zi- 



1+1* 



Z\+\* 

where we define x, y, u, v as follows: 



det 

det W 



1+1* 



i 

2 2 2 

U — V 



x 2 — y 2 



(19) 



A=A S 



X = 1 — 

y 



T 2 -fx 2 8TV(a + m 



W 

a + m) 2 



W 2 



W 



8TV 
W 



T 2 + /i 2 



u= 1 - 



W ' 



(a + m) 



W =[(a + mf + T 2 — /i 2 ] 2 + 4/i 2 T 2 



(20) 
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FIG. 1: The contours of the average phase factor R for m = 0.07. The first order phase transition 
line in 1 + 1 theory and the critical point are also shown (see discussion in Section IIIIFh . 

where a is a solution of the saddle-point equation (jl5l) . and the global minimum of Q(a) in 
Eq. (HU). 

Using Eqs. f fT9|) and fl20l) we can now obtain the average phase factor R for any values 
of T, fi and m. In general, this has to be done numerically, but in certain limiting cases, 
discussed below, explicit analytical results can be derived as well. As an illustration, the 
average phase factor contours on the T/z plane for m = 0.07 are plotted in Figure [U 



B. Arbitrary T and /i in the chiral limit m = 

In the chiral limit m = the solution of the saddle-point equation ( |T5l) can be found 
explicitly. In the high-temperature phase it is simply a = 0. Then, from ( TT9l and ( l20il the 
average phase factor R can be written explicitly: 

R-(e™) _ i(T 2 + ^) 2 -(T> + ^r 

R={e [(T2 + ^ )2 _ (T 2_ /i2)]2 ("*-«)■ ( 21 ) 

From (I2T|) . the contour where the average phase factor vanishes is determined by: T 2 +/i 2 = 1. 
The contour plot of R(T, fi) at m = is shown in Figure [2j 
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FIG. 2: The contours of the average phase factor R in the chiral limit m = 0. The chiral symmetry 
transition line and tri-critical point (TCP) are also shown. 

C. Small m and fi 



At small m and /x the solution to the fifth-order polynomial equation for the saddle point 
a can be found in the form of an expansion, which at zero temperature is given by 

(T = 0, m ~ fJ 2 < 1) , 



1 _ m fjr 
2 2 



and thus 



1 - 



2fi 2 



m 



1 - 



mi 



(T = 0, m ~ /i 2 < 1) 



(22) 



(23) 



where we used the fact that the R — value is achieved at the boundary of the pion 
condensation phase in 1 + 1* theory, and that, in QCD, the pion condensation occurs at 
fi u = —fid = to define the pion mass: m\ = m/2 in the units employed in the random 

matrix model. This result is in agreement with the earlier T = calculation [l?], fis| . 

Now, using our general result, we can extend the result of Ref. [l?], Q] to nonzero tem- 



peratures. In this case, the saddle point is given by: 



, m fl -2T 2 

a « Vl-T 2 - 

2 V 1-T 2 



4T 2 



and 



m 



2/i 2 



2 Vv 7 ! 3 ^ 2 



4/i 2 



(T < 1, m ~ /i 2 < 1) 



(24) 



1 - Vl^T 2 -q- 



(T < 1, m ~ n 2 < 1) , (25) 
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which generalizes Eq. (1231) . Equation (1251) shows that the sign problem diminishes at higher 
temperatures, which can be seen also on the contour plot in Figure [D 

D. R = contour 

The contour where R vanishes is of particular interest to our study - the sign problem 
reaches its maximum there, i.e., the fluctuations of the phase completely wash out the 
magnitude. 

We can obtain an explicit equation for the R = contour by setting u — v in equa- 
tions (USD, (120]): 

u = v [(a + to) 2 + T 2 - /x 2 ]m = 2/x 2 a . (26) 

Solving the quadratic equation ( 1261) for a and substituting the solution into the saddle-point 
equation (fl5l) one finds for the R = contour: 

T 2 1 2 ^ m2 (27) 

^ /i 2 — to 2 4(/i 2 — to 2 ) 2 ' 

in agreement with an earlier result [251 ] . 

It is interesting to consider the limiting behavior (shape) of this contour as to — > 0. As 
can be seen from Figures [Hand [2] the contour develops a singularity (a kink) at /x = 0, T = 1 
in this limit. Near this kink the shape and location of the contour scale with to, i.e., the 
contour at different values of to can be obtained by rescaling 

t -> A 2/3 t and /i 2 -> A 2/3 /i 2 as to -> Ato , (28) 

where we introduced t = T 2 — 1. This can be seen upon expanding Eq. (1271) in /x, t and m 
in the regime t : /i 2 : to 2 / 3 fixed as m — > 0: 

t = V - ^ + 0(to 4 / 3 ) . (29) 

For example, the point where the i? = contour reaches maximum temperature slides to 
T — > 1 and /x — > as 

T 2 = 1 - ^ m 2 / 3 + 0(m 4 / 3 ) and ,x 2 = ^ to 2 / 3 + 0(m^) . (30) 

This result may be useful for analysis of QCD simulations using the improved reweight- 
ing 8j or techniques which treat the sign problem by separating the phase from the absolute 



10 



value of the determinant. The reweighting methods break down at R = 0, and the problem is 
to distinguish the signatures of the critical point, which are similar to those of the breakdown 



of the reweighting 26|, [271]. By doing simulations at different values of m and comparing 
to scaling ( 1281) one can determine whether the observed signatures are those of the genuine 
critical point of 1 + 1 theory or those of the breakdown of the reweighting method. 

To be precise, in QCD, the scaling behavior of the phase transition boundary in the phase 
quenched, 1 + 1* theory (i.e., pion condensation at finite isospin chemical potential) should 
be similar to (1281) , but with critical exponent l/{(35) replacing the mean- field exponent 2/3. 
The value l/((35) ~ 0.54 is the ratio of critical scaling dimensions of the energy-like and 
ordering-field-like operators in the 0(4) universality class of the QCD phase transition (for 
two flavors). The scaling shape of the pion condensation boundary, given by Eq. ( 129]) in the 
random matrix model, in QCD will also be correspondingly different. 

E. R = domain and the sign problem 

One can understand the underlying reason that the sign problem becomes severe in the 
R = region by looking at the distribution of the zeros of the det D as a function of \x. 
These zeros can be also viewed as eigenvalues of the random matrix P, defined as 2 



(lil 2 N ~ P)7o , where 7 = • (31) 




The locations of the zeros are random, fluctuating together with matrix X in the ensemble 
in Eq. (T4J). For N — > oo the density of zeros develops finite region of support with a sharp 
boundary, which one can see already quite clearly for a finite, but large, matrix on Figure [3j 3 
The fluctuations of the phase of det D become large when [i enters the domain of the 
zeros. More explicitly, one can write 

argdetD = ^arg(/i-Ai), (32) 

i 

where the sum is over all eigenvalues A, of the matrix P. When /i is away from the domain of 
support of the eigenvalue density, the fluctuations in the eigenvalue positions do not affect 



2 This is the analog of the propagator matrix introduced in. [26 

3 At T = these distributions have been studied in Ref. [29| . 
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FIG. 3: The distribution of zeros of det D in the complex /i plane for a random matrix model with 
finite N = 1000 at T = 0.5 and m = 0.07. The filled circles indicate solutions of the R = contour 
equation (|27p . bounding the region of the maximal sign problem. 

the phase of the determinant significantly. On the other hand, when /x is in the domain 
filled with eigenvalues, there are eigenvalues which come very close to /i (as close as 1/yN) 
and even small fluctuations in the eigenvalue positions translate into large fluctuations of 
the phase, causing severe sign problem. 

To confirm this picture, one can compare the solutions of the equation fl271) with the 
distribution of the eigenvalues at the same values of T and m, as it is shown in Figure O 
The comparison clearly shows that R = (the sign problem is most severe) for the values 
of /i which fall inside the domain of eigenvalues. 



F. QCD sign problem vs RMT 

When drawing conclusions from our random matrix study for QCD one should bear in 
mind the following. In QCD, in thermodynamic V — *■ oo limit, the partition functions -Zi+i 
and Zi+i* are already different at the exponential level. Indeed, expressing the partition 
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function via the pressure, Z = exp(VP(T, //)), one finds that 4 

R = 4^ = exp[V(Pi+i - A+iO] • (33) 

As an illustration, consider sufficiently low temperatures (T <C Tn n ). With exponential pre- 
cision, the //-dependence of pressure is given by the masses of the lightest particles with 
nonzero charge to which the chemical potential couples, i.e., Pi + i(T, jj) — P 1+1 (T, 0) ~ 
/i 2 e -m N /(3T) and p 1+1 .(T,fji) - P 1+1 *(T,0) ~ fi 2 e~ m ^^ 2T \ as long as T « m n . Since 
Pi + i(T, 0) = Pi+i*(T, 0), and neglecting e~ mN ^ 3T ^ compared to e~ m7T ^ 2T ^ we can write (with 
double exponential precision): 

R(T, a*) = ~ expf-^-^] , (34) 

which means for T <C m„- the phase factor R is exponentially small for large V. 5 This is 
not surprising, since the "warm" gas of pions is very much different from the gas of baryons 
at the same temperature. 6 

In contrast, in the random matrix model, as we have seen, the exponential of the volume 
(i.e., the matrix size N) cancels in the ratio R. 

However, in QCD, for T ~ T c and higher, the difference between the pressures, 
Pi + i{T,ji) — P 1+1 *(T, n), although remaining of order Vfj, 2 , becomes smaller. Lattice studies 
indicate remarkable similarity of the partition functions Z± + i and Zi + \*, manifested, e.g., 
in the similar slopes dT c / dfi 2 of the pseudo-critical lines, and it has been also arg ued that 



the difference between the slopes is suppressed in the large iV c limit of QCD [3JJ. It may 
be also added that in QCD at asymptotically large T 3> T c the difference vanishes due 
to the asymptotic freedom - quark flavors decouple from each other, and the pressure is 
independent of the relative sign of the quark chemical potentials. Thus the random matrix 
model results might be a useful guide to the sign problem in QCD at least in the range of 
temperatures near T c , which is of much experimental and theoretical interest. 

We wish to stress again that for the comparison of QCD sign problem to the RMT to 
be meaningful the lattice 4-volume V should be finite. More quantitatively, the system 



4 For simplicity, we denote by V the 4-dimensional Euclidean volume of the finite-temperature system: 

V = F S p aC e/T. 

5 Interestingly using the condition R < 1 one could derive the mass inequality: < mjv/3 [301 ] . 

6 Of course, as T — > 0, the ratio R — > 1, which is the reflection of the fact that for T = there is no 
dependence on any chemical potential (for /i < to^/2), neither baryon nor isospin. 
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should be within the so-called epsilon-regime [32j, |33|, where the RMT description of the 
QCD becomes exact. In QCD, as the volume is increased, the role of the exponential 
factor eventually becomes dominant and R vanishes exponentially. The crossover between 
the epsilon-regime and the thermodynamic limit can be studied, e.g., along the lines of 
Ref. [34]. 

With the preceding discussion in mind, let us take the point of view [26J that the sign 
problem becomes intractable (even on a finite volume), when we enter the domain of pion 
condensation in the phase-quenched theory, i.e., R = (e 2td ) = domain of the RMT. In this 
context, it is interesting to see where, relative to this domain, is the critical point of the 
1 + 1 theory. In the 2-flavor random matrix model we studied, the critical point, as well as 
the whole first order transition line, is always inside the R = domain, as Figures Q] and [2] 
demonstrate. 



IV. SUMMARY AND CONCLUSIONS 

We investigated the strength of the QCD sign problem and its dependence on temperature 
T and quark mass m using the random matrix model. We observed that the sign problem 



diminishes at higher temperatures, which is a welcomed property, anticipated, e.g., in 20], 
and relied upon in improved reweighting techniques, e.g., Ref.[8[]. 

We also observed that the strength of the sign problem is related to the position of the 
pion-condensation region in the phase-quenched theory, equivalent to the theory where the 
baryon chemical potential is replaced by the isospin chemical potential (of the same abso- 
lute value per quark), as already discussed in Ref. 2(|. In particular, this underscores the 
importance of understanding the phase diagram of QCD at finite isospin chemical poten- 



tial W 



21 



35, 



36 



37] 



We observed and generalized to QCD the scaling behavior of the shape of the pion 
condensation region in the T\x plane near T = T c as m — > 0, where it develops a singularity 
(kink). This allows us to understand how the phase diagram of QCD with finite isospin 
chemical potential evolves towards the chiral limit. One of the practical applications of this 
result is to guard the reweighting techniques against a possible breakdown of the reweighting 
method by studying the dependence of the breakdown point on the quark mass. 

In the random matrix model we studied, the critical point of the 1 + 1 flavor theory falls 
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within the R = domain of the maximal sign problem (see Figures [T] and [2]) . This means 
that reweighting methods cannot access the critical point in such a theory. However, it is 
possible, that a 3-flavor theory, like QCD with a strange quark, where the critical point lies 
at smaller values of \x (tunable by the strange quark mass), that the critical point is outside 
the pion condensation domain. What happens in this case is an interesting question which 
we hope to address in future work. 
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